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I .  Introduction 

Compact  packaging  of  electronic  circuits  has  had  an  extraordinary  impact 
on  the  electronics  industry.  Small  scale,  integrated  circuits  allow  the 
placement  of  a  large  number  of  discrete  components  on  a  single  chip. 

However  these  small  packages  must  be  interconnected  in  order  to  perform 
specific  tasks. 

The  common  medium  for  interconnecting  these  packages  is  the  printed 
circuit  board  (PCB) .  PCB's  consist  of  a  thin,  substrate  material  on  which 
rectangular  conductors  (lands)  are  placed.  The  most  common  type  of  substrate 
or  board  is  glass-epoxy  having  a  relative  permittivity  of  approximately  e^=5. 
Typical  board  thicknesses  are  on  the  order  of  25-100  mils  with  the  typical 
thickness  being  1/16  inch  or  62.5  mils.  The  conductors  are  "etched"  on  either 
one  side  (single-sided  boards),  or  on  both  sides  (double-sided  boards)  and 
connected  through  holes  known  as  "vias".  Another  type  of  board  which  is 
somewhat  less  common  has  lands  on  one  side  and  a  ground  plane  on  the  other 
and  is  often  referred  to  as  microstrip.  There  is  also  an  increasing  use  of 
multilayer  boards  although  these  will  not  be  considered  in  this  report. 

The  metal  which  is  deposited  on  the  board  surface  is  typically  copper. 

The  thicknesses  are  characterized  by  weight  per  unit  surface  area.  The 
common  thicknesses  are  1  ounce  and  2  ounce  which  means  a  weight  of  1  or 
2  ounces  of  copper  per  square  foot.  The  1  ounce  copper  thickness  is  0.356  mm 
(1.4  mils).  The  thickness  of  2  ounce  copper  is  double  that  of  1  ounce  copper 
or  .0711  mm  (2.8  mils). 

Land  widths  vary  depending  on  the  intended  use  of  that  land.  Signal 
lands  have  widths  ranging  from  5  to  25  mils.  Lands  intended  for  carrying 
large  currents  such  as  motor  drivers  are  typcially  50  mils  in  width  and 


larger.  Land  spacings  range  from  5  mils  Co  well  over  100  mils  depending 
on  layout  constraints.  In  short,  one  can  say  very  little  about  "typical” 
land  widths  and  spacings.  In  fact,  on  most  PCB's  there  exist  areas  in 
which  no  definable  land  lengths  exist;  sections  are  filled  with  metallic 
areas  to  reduce  ground  noise  due  to  return  path  inductance.  Typical 
practice  uses  a  "gridded  ground"  to  aid  in  minimizing  "ground  shifts"  due 
to  inductive  surges  caused  by  switching  currents. 

This  wide  variability  in  land  widths,  spacing  and  shape  is  only  one 
reason  for  the  difficulty  in  characterizing  electrical  properties  of  a 
PCB.  Perhaps  the  major  problem  is  the  enormous  number  of  lands  (and  inter¬ 
connected  components)  on  a  typical  PCB.  If  one  wishes  to  model  a  PCB  for 
the  purposes  of  analyzing  functional,  electrical  performance,  one  is  faced 
with  a  modeling  as  well  as  computational  task  of  considerable  magnitude. 

Of  perhaps  equal  importance  is  the  analysis  of  the  electrical 
performance  of  a  PCB  layout  for  the  purpose  of  determining  nonfunctional 
performance  such  as  crosstalk,  ringing,  etc.  which  may  affect  functional 
performance.  Two  boards  intended  to  perform  the  same  task  but  which  have 
different  layouts  may  perform  quite  differently.  One  of  the  major  reasons 
for  this  is  crosstalk  or  electromagnetic  coupling  between  lands.  The 
current  and  voltage  associated  with  a  pair  of  lands  generate  electric  and 
magnetic  fields  which  interact  with  neighboring  lands.  These  electro¬ 
magnetic  fields  induce  signals  in  other  landr  which  appear  across  the 
terminals  of  components  attached  to  these  lands.  These  induced  signals 
may  cause  the  devices  to  malfunction  or  may  reduce  their  performance  to 
marginal  levels.  In  other  words,  on-board  interference  may  be  caused  by 
this  crosstalk. 


There  exist  a  large  number  of  digital  computer  programs  which  are 
intended  to  model  the  electrical  performance  of  interconnected  electronic 
components.  In  the  majority  of  these  codes,  the  effect  of  the  inter¬ 
connecting  lands  are  neglected.  There  is  the  inherent  ability  to  consider 
these  effects  via  lumped  circuit  approximate  models  of  the  lands  but  this 
is  rarely  done  by  the  user. 

The  speeds  of  clock  circuits  (and  the  associated  rise/fall  times)  on 
PCB's  are  increasing  at  a  dramatic  rate.  Present  clock  speeds  are  on  the 
order  of  10  to  25  mHz.  However,  speeds  of  100  mHz  are  not  far  away.  As 
the  basic  speeds  of  the  clock  circuits  increase,  so  do  the  associated 
rise/fall  times.  It  is  not  uncommon  to  find  rise/fall  times  of  pulses 
on  the  order  of  1  to  10  ns.  Subnanosecond  rise/fall  times  are  certainly 
in  the  near  future. 


As  the  speeds  and  associated  rise/fall  times  increase  so  do  their 
high  frequency  spectral  content.  For  example,  consider  a  typical  clock 
signal  approximated  as  a  periodic,  trapezoidal  waveform  shown  in  Fig.  l-l(a) 
The  repetition  frequency  is  denoted  as  f  with  the  bit  rate  or  period  being 
T  =  1/f  .  The  pulse  width  (between  50%  points)  is  denoted  as  t.  The 
rise  and  fall  times  are  denoted  as  and  t ^ ,  respectively.  For  equal 
rise  and  fall  itnes,  one  can  obtain  the  Fourier  series  coefficients  as 


.  /  n7IT  \  .  ,n  71  "V  \ 
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o  o 


(1.1) 


The  amplitudes  of  these  spectral  components  can  be  bounded  as  shown  in 
Fig.  1.1(b).  The  spectra  are  approximately  constant  out  to  a  frequency 


of  l/rt.  Between  this  and  1/ttt^  they  fall  off  at  -20  dB/decade.  Above 


3 


l/nx  they  fall  off  at  -40  dB/decade.  Thus  the  high-frequency  spectral 
content  of  the  signal  is  primarily  determined  by  the  rise/fall  time  of  the 
pulse . 

For  example,  consider  a  typical  10  mHz  signal  with  50%  duty  cycle  having 
5  ns  rise/fall  times.  The  breakpoints  are  at  1.6  MHz  and  63.7  MHz.  A  pair 
of  lands  of  typical  length  of  10  cm  on  a  glass  epoxy  board  will  be  1/10  of  a 
wavelength  and  electrically  short  at  approximately  170  mHz.  Thus  the  lands 
will  be  electrically  short  at  frequencies  below  the  breakpoints.  In 
this  range,  the  crosstalk  typically  increases  at  +20  dB/decade.  The  spectrum 
of  the  crosstalk  signal  can  then  be  found  as  the  product  of  the  spectrum  of 
the  clock  signal  and  the  spectrum  of  the  crosstalk  transfer  function  as  shown 
in  Fig.  1.2.  Note  that  the  resulting  spectrum  of  the  crosstalk  pulse  is 
constanf  between  the  above  two  frequencies.  Reducing  the  pulse  rise/fall 
time  directly  increases  the  high  frequency  content  of  the  crosstalk  spectrum. 
It  will  be  shown  that  this  results  in  an  increased  amplitude  of  the  time- 
domain  crosstalk  pulse.  Thus  the  potential  interference  of  this  crosstalk 
pulse  increases  directly  as  the  rise/fall  time  of  the  source  pulse  decreases. 

This  concept  will  be  addressed  in  more  detail  in  later  chapters. 

However  it  is  sufficient  to  point  out  here  that  potential  crosstalk  problems 
will  increase  with  anticipated  pulse  rise/fall  time  reductions.  Crosstalk 
which,  with  present  speeds,  cause  few  problems  will  generate  an  ever 
increasing  set  of  interference  problems  in  the  future. 

It  is  therefore  becoming  increasingly  important  to  consider  crosstalk 
on  PCB's.  Present  computer-aided  design  (CAD)  programs  do  not  generally 
include  modeling  capabilities  for  the  purpose  of  including  crosstalk. 

There  are  inherent  capabilities  via  lumped  circuit  approximations  of  the 


land-circuits  as  will  be  discussed  in  later  chapters.  However,  a  circuit 
designer  who  uses  these  codes  to  check  functional  performance  of  a  board 
will  not  be  predisposed  to  include  crosstalk  considerations  because  pre¬ 
paration  of  crosstalk  models  for  use  in  the  program  will  require  a  large 
effort  on  the  analyst's  part.  For  example,  the  analyst  will  generally 
not  have  information  readily  available  as  to  the  land  dimensions  or 
configuration;  functional  schematics  do  not  include  this.  Thus  the  first 
major  problem  which  the  analyst  will  face  if  he/she  intends  ro  include 
a  crosstalk  analysis  in  the  functional  performance  analysis  is  gathering 
these  physical  layout  data.  Generally  these  data  are  computerized  since 
the  board  layout  was  probably  done  by  automated  design  techniques. 

Once  the  analyst  gathers  this  physical  layout  data,  he/she  is  faced 
with  a  more  formidable  problem;  translating  that  physical  data  into 
mathematical  models  suitable  for  use  in  the  intended  CAD  program.  If 
gathering  the  physical  layout  data  is  not  a  sufficient  deterent  to  including 
crosstalk  analysis  in  the  CAD  program,  this  latter  aspect  of  the  problem — 
preparing  a  mathematical  model — generally  is.  What  would  be  highly  desir¬ 
able  if  we  desire  crosstalk  analysis  to  be  integrated  into  functional 
analysis  would  be  to  make  the  inclusion  of  these  crosstalk  models  as 
simple  as  possible  for  the  analyst. 

The  primary  purposes  of  this  report  are  (1)  to  examine  the  available 
time-domain,  crosstalk  prediction  models  and  to  determine  their  suitability 
for  inclusion  in  CAD  programs  and  (2)  to  develop  a  model  which  minimizes 
the  implementation  difficulty  while  providing  reasonable  estimates  of 
FOB  crosstalk  and  is  suitable  for  inclusion  in  typical  CAD  programs.  It 
is  certainly  not  possible  to  determine  an  all-inclusive  prediction  model 


which  is  simple  and  at  the  same  time  vl.  handle  general,  transmission 
lines  of  any  dimensions  driven  by  i  ei  e  lv  general  signals.  We  will 
evaluate  models  using  criteria  such  as  tvpical  board/land  dimensions  and 
typical  driving  signals.  Within  these  limitations  we  will  find  that  one 
can  develop  fairly  simple  models  which  yield  quite  accurate  predictions. 
Within  these  restrictions  we  will  find  that  the  only  limitation  to 
realizing  our  goal  of  making  crosstalk  analysis  a  simple  matter  in  CAD 
programs  is  translating  the  physical  layout  dimensions  to  the  appropriate 
per-unit-length  line  parameters. 

1.1.  The  Conventional,  Transmission-Line  Model 

The  usual  method  for  modelling  electromagnetic  coupling  between  parallel 
conductors  is  the  distributed-parameter,  transmission  line  model  [1].  The 
model  assumes  that  the  mode  of  propagation  on  the  line  is  the  Transverse 
Electro  Magnetic  (TEM)  mode.  For  lines  which  have  (1)  imperfect  conductors 
and/or  (2)  an  inhomogeneous  surrounding  medium,  a  pure  TEM  mode  is  not 
possible.  However  the  dominant  mode  of  propagation  resembles  the  TEM  mode 
and  is  referred  to  as  the  "quasi-TEM"  mode.  Assuming  a  TEM  mode  of 
propagation  for  these  structures,  while  technically  not  correct,  never¬ 
theless  has  been  shown  to  provide  accurate  representation  for  (1)  typical 
conductor  materials  and  (2)  typical  dielectric  inhomogeneities  so  long 
as  the  frequency(s)  of  excitation  is  small  enough.  For  typical  PCB's, 
the  transmission  line  model  will  provide  an  accurate  characterization 
for  frequencies  up  to  the  GHz  range.  Therefore  the  fundamental  model  which 


we  will  use  is  the  distributed  parameter  transmission-line  model. 


Consider  a  general  transmission  line  consisting  of  (n+l)  lossless 

conductors  immersed  in  a  lossless,  surrounding  medium.  We  will  assume 

the  line  to  be  uniform;  that  is,  cross-sectional  views  at  any  two  points 

along  the  line  will  be  identical  in  dimensions  as  well  as  conductor  and 

medium  properties.  Thus  we  consider  (n+l)  parallel  conductors.  The 

conductors  will  be  parallel  to  the  x  axis.  An  electrically  small  Ax 

length  of  the  line  is  modeled  as  shown  in  Fig.  1.3.  with  per-unit-length 

self  inductances  (capacitances)  1..  (c..)  and  mutual  inductances 

n  11 

(capacitances)  (c„)  ^or  i « j= * » • • * »n •  These  per-unit-length  parameters 

are  calculated  from  static  (DC)  configurations  of  the  cross-sectional 
fields  [1],  Line  voltages,  V.  (x,t),  are  defined  for  n  of  the  conductors 
with  respect  to  the  (n+l)st  (reference)  conductor.  We  refer  to  this  con¬ 
ductor  as  the  0  conductor.  Line  currents , ^^(x,t)  are  defined  for  the  n 
conductors  in  the  positive  x  direction  and  return  in  the  reference  conductor 
so  that  the  net  current  in  the  +x  direction  at  each  point  along  the  line 
is  zero.  These  notions  are  exact  for  the  TEM  mode  and  are  reasonable 
approximations  for  the  "quasi-TEM"  mode  [1]. 

From  the  circuit  in  Fig.  1.3  one  can  obtain  the  transmission-line 
equations,  in  the  limit  as  Ax  ->  0  as  a  coupled  set  of  2n  partial  differ¬ 
ential  equations  [1]: 

aiAx.t)  a 

— 5 -  =  -  L  ^S(x,t)  (1.2a) 

dx  ~  dt  — 

9^(  x ,  t )  dV(  x ,  t ) 

- 5 -  =  -  C  — 5- -  (1.2b) 

dx  ~  dt 
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*nAX 


Jl  j (x+AX,t  ) 


Reference  Conductor 


where  the  complex,  phasor  voltages  and  currents  are  denoted  as  V^(x)  and 
Ii(x),  respectively,  and  to  =  2irf  is  the  radian  frequency  of  excitation. 
For  this  case,  the  transmission  line  equations  become  [1] 

V(x)  -  -jtoL  I_(x)  (1.6a) 

I(x)  =  - jtoC  V(x)  (1.6b) 

where  [V(x)].  =  V.(x)  and  [ X ( x ) ] .  =  I.(x). 

—  i  l  —  i  i 

For  a  homogeneous  surrounding  medium  (which  is  not  the  case  for 
PCB's)[l] 


L  C  =  uel 
~  ~  ~n 


(1.7) 


where  1^  is  the  nxn  identity  matrix  and  y(e)  is  the  permeability 
(permittivity)  of  the  surrounding  medium.  The  velocity  of  propagation 
in  this  homogeneous  medium  is 


(1.8) 


/\ie 

Equation  (1.7)  shows  that  we  need  only  determine  one  parameter  matrix.  Fo 
example,  if  we  obtain  the  per-unit-length  capacitance  matrix  C  from  static 
field  considerations,  then  the  per-unit-length  inductance  matrix  is  easily 
obtainable  in  terms  of  the  inverse  of  the  capacitance  matrix  as 


L 


(1.9) 


v 

For  an  inhomogeneous  medium  we  can  obtain  both  parameter  matrices 


from  only  capacitance  calculations  as 


where  is  the  per-unit-length  capacitance  matrix  with  the  dielectric 

insulation  removed  and  v  =  l//p  £  =  3x10®  m/s  is  the  velocity  of 

o  o  o 

propagation  in  free  space. 

The  matrices  L  and  C  contain  all  of  the  information  concerning  the 
cross-sectional,  physical  configuration  of  the  line.  Tor  the  purposes  of 
this  report,  we  will  assume  that  these  parameters  have  been  calculated  (or 
measured)  and  are  known.  The  primary  thrust  of  this  report  is  to  develop 
a  simple  solution  of  (1.2)  suitable  for  implementation  in  CAD  programs. 


1.1.2.  The  Transmission-Line  Equations  -  Lossy  Lines 

The  previous  formulation  assumed  perfect  conductors  and  a  lossless 
medium.  Losses  in  the  surrounding,  inhomogeneous  medium  pose  no  problems 
in  so  far  as  representation  in  the  per-unit-length  model  of  Fig.  1.3.  In 
order  to  consider  a  lossy  surrounding  medium  we  add  per-unit-length 
conductances  between  each  of  the  n  conductors  and  between  each  of  the  n 
conductors  and  the  reference  conductor.  Similarly,  imperfect  conductors 
may  be  considered  in  the  per-unit-length  model  by  adding  per-unit-length 
impedances  in  series  with  each  conductor.  For  typical  PCB's  losses  in 
the  medium  are  neglectable  for  frequencies  below  1  GHz.  However,  conductor 
losses  may  be  significant  for  all  frequencies.  A  simple  example  is  the 
resistance  of  the  reference  conductor.  Currents  of  the  other  circuits 
pass  through  this  reference  conductor  and  the  resulting  voltage  drop 
across  this  reference  conductor  appears  directly  in  all  the  other  circuits. 
This  phenomenon  is  referred  to  as  "common  impedance  coupling"  or  "ground 
loop  coupling"  and  may  be  important  at  very  low  frequencies.  In  ribbon 
cables  where  all  circuits  share  a  common  return  wire,  this  component  of 


crosstalk  coupling  has  been  shown  to  be  dominant  below  frequencies  as  high 
as  10  kHz  [2,3].  Thus  conductor  losses  are  not  necessarily  neglectable. 

However,  if  we  include  imperfect  conductors,  we  are  faced  with  several 
formidable  modeling  problems.  The  first  of  these  is  determining  a  suitable 
model  for  characterizing  these  losses.  In  the  case  of  sinusoidal,  steady- 
state  excitation  of  the  line,  it  is  a  rather  straightforward  task  to  model 
the  effects  of  conductor  losses;  simply  add  frequency  dependent  impedances 
in  series  with  each  conductor  [ 1 j .  These  impedances  are  frequency  dependent 
as  a  result  of  skin  effect  and  vary  typically  as  the  square  root  of  frequency 
Thus  we  add 


Vf)  -  Zi0  *zu/I 


(1.11) 


in  series  with  each  conductor  (including  the  reference  conductor).  The 
constants  and  Z^  are  dependent  on  the  conductor  cross-section  and 
physical  properties.  Similarly  losses  in  the  medium  are  frequency  dependent 
and  may  be  modeled  by  adding  a  frequency  dependent  conductance,  g^(f) 
between  each  conductor.  This  results  in  modifying  the  transmission-line 
equations  for  sinusoidal  steady-state  excitation  in  (1.6)  as 


where 


V(x)  -  -CZ( f )  +  jwL]  I ( x ) 


I(x)  -  — [G(  f )  ♦  jcuC]  V(x) 


[Z( f)3 . .  -  Z.(f)  +  Zn(f) 

~  li  i  0 


[Z(f)]..  -  z0(f) 


(1.12a) 


(1.12b) 


(1.13a) 


(1.13b) 
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where  Zg(f)  is  Che  per-unit-lengCh  impedance  of  the  reference  conductor 


(1.13) 


Including  losses  in  the  transmission-line  model  is  therefore  not 
conceptually  difficult  in  the  case  of  sinusoidal,  steady-state  excitation. 

Of  course,  determining  the  constants  Z^q  and  Z^j  in  (1.11)  and  the  function 
g^j(f)  in  (1.13)  requires  analysis  but  once  that  is  concluded,  there  is 
no  difficulty  in  including  these  effects  in  the  transmission-line  model. 

On  the  other  hand,  including  lossy  conductors  and  media  in  the  time- 
domain,  transmission-line  model  and  the  resulting  equations  in  (1.2)  presents 
considerable  difficulty.  For  example,  how  does  one  represent  frequency 
dependent  losses  Z^(f)  and  g^(f)  in  the  time  domain?  The  usual  method 
of  doing  this  is  to  omit  the  frequency  dependence  of  the  losses  and  assume 
constant,  per-unit-length  resistances  in  series  with  each  conductor  and 
conductances  between  each  line.  The  time-domain,  transmission-line 


equations  in  (1.2)  become  [1] 


3l/(x,t) 


-  [R  +  jU)L] 


^A(x,t) 


where 


3$(x,t)  3\i(x,t) 
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(1.14b) 


(1.15a) 


(1.15b) 


(1.15c) 
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and  these  parameters  are  constants. 


Solution  of  the  time-domain  transmission-line  equations  without  losses 
has  been  considered  in  numerous  instances  [4-8].  Inclusion  of  frequency 
independent  losses  as  in  (1.14)  has  been  the  usual  method  of  including 
losses  in  time-domain  solutions  [9-11].  However  frequency  independent 
losses  as  in  (1.14)  are  a  gross  simplification.  Frequency  independent 
losses  are  usually  assumed  simply  because  they  are  easier  to  model  than 
are  frequency  dependent  losses.  A  recent  attempt  at  including  frequency 
dependent  losses  in  the  time-domain  model  is  given  in  [12],  However  this 
requires  measurements  of  the  line  under  study  and  as  such  is  not  satisfactory 
as  a  prediction  model. 

A  recent,  promising  approach  to  modeling  frequency  dependent  losses 
in  the  time-domain  model  is  given  in  [13].  This  model  uses  a  lumped 
circuit  to  approximate  the  /f  impedance  dependence  of  the  conductors  as 
shown  in  Fig.  1.4.  The  values  of  the  resistors  and  inductors  are  chosen 
such  that  as  frequency  is  increased  each  branch  "opens  up"  thereby  increasing 
the  resistance  in  such  a  way  that  it  becomes  a  function  of  /f  up  to  an 
upper  frequency  limit.  For  example,  consider  m=3  in  Fig.  1.4.  The  input 
admittance  is 


Y 


ju)L2  + 


2  1 

R,  +  — 


3  ju>L, 


(1.16) 


At  DC  the  impedance  is  the  parallel  combination  of  ,  R2>  R^.  As  frequency 
is  increased  the  limiting  value  is  R^ .  For  intermediate  frequencies  the 
impedance  approximates  a  /f  dependence  [13]. 


The  advantage  of  this  circuit  is  that  it  can  be  included  in  a  time- 


domain  CAD  model  directly  with  no  translation  between  frequency  domain 
dependence  and  time-domain  dependence.  Nevertheless  this  model  does  not 
show  how  the  time-domain,  transmission-line  equations  in  (1.2)  can  be 
modified  to  include  frequency  dependent  losses.  Thus  inclusion  of  frequency 
dependent  losses  in  the  time-domain,  transmission-line  equations  is  still 
an  open  problem. 

1.2.  Incorporating  the  Terminal  Constraints 

The  exact  solution  to  the  transmission  line  equations  in  (1.2)  is 
(even  without  regard  to  implementation  in  CAD  programs)  not  a  simple  matter. 
An  exact  solution  will  be  shown  in  the  next  chapter.  Traditionally  solutions 
to  these  equations  are  obtained  via  approximations  also  discussed  in  the 
next  chapter. 

Nevertheless  the  general  solution  to  the  transmission  line  equations 
involve  4n  undetermined  constants.  The  final  step  in  the  complete  solution 
to  these  equations  is  the  incorporation  of  the  terminal  constraints  at 
the  two  ends  of  the  line.  Consider  a  line  of  length  9£  extending  from  x*0 
to  x-/.  At  x*0,  the  voltages  and  currents  are  related  or  constrained  in 
some  fashion: 

f^(v(0  ,t) ,  1(0, t)  =  0  (1.17a) 

Similarly  at  x  -  X 

f^vQC.t),  l(X,t>)  -  0  (1.17b) 

If  the  terminations  are  linear  (as  they  generally  are  not  on  PCB's)  then 
(1.17)  can  be  represented  as  linear  constraints.  For  sinusoidal,  steady- 
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state  excitation,  linear  termination  networks  are  usually  represented  as 
generalized  Thevenin  or  Norton  equivalents  [l].  Nevertheless,  the  complete 
solution  to  the  transmission-line  equations  involves  incorporating  the 
line  terminal  constraints  into  the  general  solution.  This  final,  important 
part  of  the  solution  process  will  be  the  most  difficult  part.  Nonlinear 
termination  networks  as  is  usually  the  case  for  PCB's  complicate  the  process 
considerably  simply  because,  for  nonlinear  terminations,  we  cannot  write 
simple,  terminal  constraint  equations. 


1.3.  Report  Summary 

Chapter  II  will  be  devoted  to  a  discussion  of  previous  methods  for 
direct,  time-domain  solution.  The  advantages  as  well  as  the  limitations 
and  disadvantages  for  implementation  in  CAD  codes  will  be  discussed. 

In  Chapter  III,  a  simple,  time-domain  model  which  is  suitable  for 
inclusion  into  CAD  codes  will  be  developed.  Although  an  approximation, 
the  simplicity  of  the  model  will  allow  the  analysis  of  large,  practical 
PCB's.  The  various  limitations  and  restrictions  of  this  model  will  be 
investigated  fully. 

Experimental  re  llts  will  be  given  in  Chapter  IV  which  show  the 
prediction  accuracy  of  the  simple  model  as  well  as  its  limitations. 

A  summary  of  the  report  will  be  given  in  Chapter  V  along  with 
recommendations  for  future  investigations. 
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Survey  of  Previous  Methods 


The  purpose  of  this  chapter  is  to  summarize  the  methods  for  direct, 
time-domain  analysis  of  transmission  lines  for  the  purpose  of  predicting 
crosstalk.  Each  of  these  methods  will  be  discussed  from  the  standpoint 
of  their  limitations,  advantages  and  disadvantages  and  their  suitability  for 
inclusion  in  CAD  codes. 

It  should  be  kept  in  mind  that  the  fundamental  model  of  the  transmission 
line  which  all  models  strive  to  approximate  is  the  TEM  mode,  distributed 
parameter  model  discussed  in  the  previous  chapter  and  represented,  for 
the  lossless  case,  by  equations  (1.2).  Therefore  there  are  no  "new 
theoretical"  models  but  only  methods  for  solving  or  approximating  the 
solution  to  the  transmission-line  equations. 


2.1.  Lumped  Circuit  Iterative  Models 


Lumped  circuit  iterative  models  are  the  most  common  forms  of  models 
of  transmission  lines  [l].  Of  all  the  various  methods  which  are  used 
in  CAD  codes  for  modeling  transmission  lines,  the  lumped  circuit  iterative 
models  are  the  most  common.  There  exist  only  a  few  exceptions  to  this. 

The  philosophy  of  lumped  circuit  iterative  models  is  as  follows. 

It  is  well  known  that  lumped  circuit  elements  are  valid  representations 
of  an  object  so  long  as  the  physical  dimensions  are  electrically  small. 

The  basic  TEM  mode  model  of  a  line  assumes  that  the  cross-sectional  dimen¬ 
sions  of  the  line  are  electrically  small  otherwise  higher  order  modes  are 
possible.  Thus  the  transmission  line  is  divided  into  N  sections  each  of 
which  is  electrically  short  at  the  frequency  of  interest  as  shown  in 
Fig.  2.1.  Traditionally  the  line  is  divided  into  N  equal  length  sections 
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each  of  length  3-/N  such  that 


where 


3f/N  «x 


X  =  v/f 


(2.1) 


(2.2) 


is  a  wavelength  at  the  frequency  of  interest. 

Once  this  is  done,  each  section  is  modelled  with  a  lumped  circuit 
model  to  represent  the  distributed  parameter  model  of  the  line  in  Fig.  1.3. 
There  are  typically  four  topologies  which  one  finds.  These  are  arrived 
at  by  either  placing  all  parallel  elements  on  the  left  or  on  the  right 
of  the  series  elements  or  placing  half  the  parallel  elements  on  the  left 
and  right  of  the  series  elements  or  placing  half  the  series  elements  on  the 
left  and  right  of  the  parallel  elements  and  are  referred  to  as  the  lumped 
r,*l,  Pi  and  Tee  circuits,  respectively.  These  are  shown  in  Fig.  2.2,  2.3, 
2.4  and  2.5. 

The  various  approximation  capabilities  of  these  lumped  circuit 
iterative  models  was  investigated  in  [14].  For  lossless  lines,  one  can 
characterize  the  performance  of  these  models  in  terms  of  the  ratio  3?/X 
rather  than  frequency.  It  was  found  that  one  section  of  each  model  provided 
an  accurate  approximation  of  the  transmission  line  for  frequencies  up  to 
a  point  where  the  line  is  approximately  X  long.  Adding  sections 
increases  the  coverage  in  a  nonlinear  fashion.  For  example,  it  was  found 
that  10  sections  do  not  uniformly  increase  the  prediction  accuracy  to 
frequencies  where  the  line  is  one  X  long. 

In  addition,  it  was  found  that  the  prediction  frequency  range  was 
dependent  on  the  termination  impedances.  For  example,  suppose  the 
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Fig.  2.3.  The  lumped  *"1  circuit 


termination  impedances  are  small  in  value.  Then  the  capacitances  on 
the  right  of  the  Pi  model  are  effectively  of  no  consequence  and  the  Pi 
model  essentially  reduces  to  the  H  model. 

Thus  the  precise  frequency  range  for  which  lumped-circuit  iterative 
models  provide  an  adequate  approximation  to  the  transmission  line  model  is 
(1)  very  difficult  to  predict,  (2)  dependent  on  the  terminations,  and  (3) 
not  significantly  increased  with  the  addition  of  more  sections.  These 
deficiencies  make  the  lumped-circuit  iterative  models  difficult  to  apply 
in  CAD  codes.  In  addition,  if  the  line  is  electrically  long  at  the  fre¬ 
quency  of  interest,  a  large  lumped  circuit  is  needed  to  model  the  line. 

For  PCB's  which  contain  a  very  large  number  of  coupled  lands,  one  is 
essentially  restricted  to  model ing  frequencies  where  all  the  traces  on  the 
board  (or  most  of  them)  are  electrically  short.  In  the  case  of  pulse 
excitation,  such  as  clock  or  trapezoidal  pulses  discussed  in  Chapter  I, 
one  is  restricted  to  rise/fall  times  for  which  the  second  break  point  in 
Fig.  1.1,  1/ttt^,  is  much  less  than  the  frequency  at  which  the  lands  are, 
say,  X  in  length.  For  a  typical  PCB  land  length  of  10  cm  one  can 
consider  rise/fall  times  no  shorter  than  10  to  20  ns.  In  the  next  chapter 
we  will  consider  a  much  simpler  model  which  has  essentially  the  same 
restriction. 

However,  the  main  advantage  of  the  lumped-circuit  iterative  models 
is  that  they  require  no  understanding  of  specialized  solution  methods  and, 
given  the  per-unit-length  capacitances  and  inductances,  one  can  readily 
implement , the  models  in  a  CAD  program  using  lumped-circuit  elements  which 
are  already  available  in  the  CAD  program  in  a  fashion  very  similar  to  the 
process  the  analyst  is  presently  using  to  model  the  functional  modules. 
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Although  there  are  several  such  transformations  which  have  been  used 


[5,6]  an  example  is 


T 

~V 


L  T 


=  T 


where 


T-1  C  L  T  -  A 


(2.7a) 


(2.7b) 


(2.8) 


is  diagonal.  The  columns  of  T  are  the  eigenvectors  of  the  matrix  product 
C  L,  and  the  entries  along  the  main  diagonal  of  A  are  the  associated 
eigenvalues  [1].  For  this  transformation,  (2.5)  becomes 

J^-V^x.t)  ■■  -  J^-*$(x,t)  (2.9a) 

dx  —  dt  — 

^■^}(x,t)  ““A^-V^x.t)  (2.9b) 

Thus  we  have  reduced  n  coupled  lines  to  a  set  of  n  uncoupled,  two-conductor 
lines  each  of  which  are  characterized  by 


3 

3x 


l/(x,t) 


g^.U.t) 


(2.10a) 


Thus  we  have, 
two-conductor 


,t) 


(2.10b) 


in  terms  of  these  modal  variables,  n  sets  of  uncoupled, 
lines  each  of  which  have  a  characteristic  impedance  of 


(2.11) 
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and  velocity  of  propagation 


Vi  *  nr 

i 


(2.12) 


as  shown  in  Fig.  2.6. 


Consider  a  two-conductor,  lossless  line  with  characteristic  impedance 


Z  and  velocity  of  propagation  v.  The  transmission-line  equations  are 

L. 


3l/(x,t)  „  3^(x,t) 

3x  "  3 


3^(x,t)  3l/(  x ,  t ) 

3x  *  c  3t 


where 


N 


/  Zc 


The  general  solution  to  these  equations  is 


l/(x,t)  -  l/(t-x/v)  +  l/(t+x/v) 


^(x,t)  =•  V+(t-x/v)  -  yf  (t+x/v) 

ZC  zc 


At  the  terminals  we  have 


1/(0, t)  =  l/(t)  +  i?(t) 


<$(0,t)  =  ^--/(t)  -  V(t) 

zc  zc 


V(£,t)  -  l/(t-T)  +  V(t  +  x) 


^S(i.t)  -  -  ^-V(t  +  T) 


(2.13a) 


(2.13b) 


(2.14a) 


(2.14b) 


(2.15a) 


(2.15b) 


(2.16a) 


(2.16b) 


(2.17a) 


(2.17b) 
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where  the  one-way  line  delay  is 


T  -  £/v  (2.18) 

Solving  (2.16)  and  (2.17)  yields 

]/(0,t)  +  Zc^(0.t)  =l/(^,t-T)  +  Z^tf.t+l)  (2.19a) 

V(0,t)  -  Zct6(0,t)  =V(*,t-T)  -  Z(?S(4t-T)  (2.19b) 

l/(£,t)  +  Z^tf.t)  -  l/(0,t-l)  +  Z^SKO.t-i)  (2.19c) 

,t)  -  Z^tf.t)  -  T/to.t+x)  -  Zc-A(0,t+T)  ( 2 . 19d) 

Equation  (2.19b)  yields 

V(0,t)  =  [l/(/,t-T)  -  Z^Of.t-T)]  +  Zc«l(0,t)  (2.20) 

Similarly  equation  (2.19c)  yields 

/(Z.t)  -  [-/(O.t-T)  +  Zc^(0,t-t)]  -  ZcJ(*,t)  (2.21) 


This  yields  the  equivalent  circuit  shown  in  Fig.  2.7  involving  controlled 
sources  with  time  delay  as  first  proposed  by  Branin  [15]  and  later  extended 
to  include  frequency  independent  losses  by  Liu  [16].  This  model  has  been 
used  by  a  host  of  others  for  two-conductor  lines  in  one  form  or  another  [17] 
and  extended  to  multiconductor  lines  with  the  above  modal  transformation  by 
Chang  [5].  The  technique  gives  an  extremely  efficient  numerical  computation 
technique  [15].  Its  principle  disadvantage  is  in  CAD  implementation  for 
multiconductor  lines  which,  according  to  the  above,  requires  determination 
of  the  modal  transformation  matrices  (eigenvectors)  of  certain  matrices  as 
in  (2.6).  Thus  implementation  of  this  technique  for  multiconductor  lines 


can  be  very  time  consuming,  computationally,  even  for  small  numbers  of 
coupled  lines. 


In  the  case  of  direct,  time-domain  analysis  of  lossy  lines,  we  are 
faced  with  the  usual  problem  of  how  to  include  these  frequency  dependent 
losses  in  the  time-domain  transmission  line  equations. 

2.3.  Frequency-Domain  Methods 

Perhaps  the  method  for  time-domain  analysis  which  is  second  in 
popularity  to  the  lumped-circuit  iterative  models  is  referred  to  as  the 
frequency  domain  method.  First  one  determines  (through  calculation  or 
measurement)  the  frequency-domain  transfer  function  (magnitude  and  phase) 
relating  the  input  signal  to  the  desired  crosstalk  signal,  T( f ) .  Then  one 
determines,  via  Fourier  methods,  the  spectrum  (magnitude  and  phase)  of 
the  input  pulse.  Multiplying  the  magnitude  functions  and  adding  the 
phase  functions  gives  the  spectrum  (magnitude  and  phase)  of  the  crosstalk 
pulse  as  shown  in  Fig.  2.8.  The  spectrum  of  this  crosstalk  pulse  is  then 
transformed  via  the  inverse  Fourier  transform  to  yield  the  time-domain 
crosstalk  pulse. 

This  method  is  quite  popular  but  not  suited  to  implementation  in  CAD 
codes.  A  major  advantage  of  this  method  is  that  frequency-dependent  losses 
are  readily  modeled.  A  major  disadvantage  is  that  the  transfer  function 
must  represent  a  linear  system  since  the  method  inherently  relies  on 
the  principle  of  superposition  to  process  the  individual  spectral  components 
of  the  input  pulse.  Thus,  although  the  line  may  be  linear,  an  important 
part  of  the  overall  transfer  .  inction  is  the  line  terminations.  If  these 
are  nonlinear,  the  overall  system  is  nonlinear  and  this  method  cannot  be 
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used.  Since  Che  terminations  on  PCB  lands  are  electronic  modules  which  are 
inherently  nonlinear  (highly  nonlinear  in  their  digital  mode  of  operation), 
this  method  is  not  a  viable  one. 

2.4.  Summary 

Of  all  the  models  considered,  only  the  lumped-circuit  iterative  models 
are  appropriate  for  use  in  CAD  codes  when  multiconductor  lines  are  involved 
A  number  of  numerical  computation  techniques  have  been  proposed  in  the 
literature  for  analyzing  "transmission  lines"  directly  in  the  time  domain 
[15-18].  However,  with  the  exception  of  the  modal  method  proposed  by  Chang 
none  of  these  methods  are  capable  of  handling  multiconductor  lines,  i.e., 
predicting  crosstalk. 

The  predictive  ability  of  the  lumped-circuit  iterative  models  will 
be  investigated  in  Chapter  IV  with  reference  to  experimental  results. 
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A  Simple,  Time-Domain  Prediction  Model 


Methods  for  the  solution  of  the  transmission-line  equations  were  given 
in  the  previous  chapter.  Only  one  method — the  lumped-circuit,  iterative 
model — was  suitable  for  inclusion  into  CAD  codes.  The  lumped-circuit, 
iterative  model  relied  on  dividing  the  line  into  electrically  short  lengths 
and  modeling  those  sections  with  lumped  circuits  which  resembled  the 
per-unit-length  model  from  which  the  transmission  line  equations  were 
derived . 

The  primary  advantage  of  this  technique  is  that,  given  the  per-unit- 

length  parameters  of  the  line,  it  is  a  simple  matter  to  incorporate  these 

models  into  lumped-circuit,  CAD  codes  using  the  available  model  elements 

of  the  code  (resistors,  inductors,  mutual  inductors  and  capacitors).  The 

primary  disadvantage  of  this  method  is  that  one  must  divide  the  line  into 

electrically  short  sections  and  model  each  section  with  a  fairly  complex 

lumped  circuit.  If  the  entire  line  is  electrically  short  (say,  less  than 
J_ 

10  A  in  length),  then  one  section  will  usually  suffice.  However,  further 
segmentation  of  the  line  usually  does  not  increase  the  frequency  coverage 
significantly.  Even  if  one  could  satisfactorily  model  a  line  by  dividing 
it  into  sections  10  A  in  length,  a  one  wavelength  section  of  line  would 
require  10  sections  which  is  a  very  complex  circuit.  Modeling  all  lands 
on  typical  PCB's  in  this  fashion  will  generally  result  in  an  overall, 
lumped  circuit  which  is  prohibitively  complex  either  in  computation  time 
or  required  array  storage.  So,  even  with  the  lumped  circuit  iterative 
models,  one  is  typically  restricted  to  electrically  short  lines. 


If  we  consider  typical  PCB  land  dimensions  and  pulse  rise/fall  times, 
it  is  found  that  requiring  the  lands  to  be  electrically  short  is  not  a 
serious  restriction.  For  example,  suppose  we  consider  trapezoidal  pulse 
signals  such  as  are  shown  in  Fig.  1.1  and  are  typical  of  clock  signals. 
Suppose  we  consider  10  ns  rise/fall  times.  The  second  break  frequency  in 
the  spectrum  occurs  at  31.8  mHz.  Suppose  we  require  that  the  spectrum  decay 
by  10  dB  at  some  point  after  this  second  break  frequency.  This  point 
becomes  100  mHz.  If  we  now  consider  land  lengths  which  are  no  longer  than 

JL 

10  X  at  100  mHz,  this  means  that  we  are  restricted  to  analyzing  land 
lengths  no  longer  than  15  cm.  If  the  pulse  rise/fall  time  is  reduced  to 
1  ns,  the  maximum  line  lengths  reduce  to  1.5  cm.  Therefore  typical  board 
dimensions  will  be  capable  of  being  analyzed  with  one  lumped-circuit 
iterative  model  per  set  of  lands  so  long  as  the  pulse  rise/fall  times  are 
not  much  less  than  10  ns. 

In  the  above  frequency-domain  criterion  we  are  assuming  that  the  high 
frequency  components  of  the  spectrum  (those  above  the  second  break  point, 

1 / TTTr )  have  decayed  sufficiently  at  some  upper  point,  f^,  such  that  incorrect 
prediction  of  them  will  result  in  negligible  errors  in  the  prediction  of 
the  crosstalk  pulse.  This  can  be  equivalently  related  to  the  time  domain. 

For  example  suppose  that  we  require  that  the  line  be  electrically 
short  at  some  frequency  f 


For  example,  choose  k  ■  yg-  .  Next,  suppose  that  this  frequency  is  well 
above  the  second  break  point  of  the  input  pulse  spectrum: 

f  >  -L_ 

u  -  ITT  0.2) 

r 


For  a  fall-off  of  10  dB  from  the  level  of  1/ttt^,  we  require  a  -  3.16.  The 
one-way  time  delay  of  the  line  is 


T 


d 


(3.3) 


Substituting  (3.3)  and  (3.2)  into  (3.1)  yields 


>  a 

T  T  -  -  T 

r*  f  kit  d 


(3.4) 


Thus  the  rise/fall  time  of  the  pulse  must  be  much  longer  than  the  one-way 
delay  time  of  the  line.  For  k  -  ,  a  *  3.16, 

T  T.  >  10  T.  (3.5) 

r,  f  -  d 

This  result  shows  that,  so  long  as  the  pulse  rise/fall  time  is,  say,  a  factor 
of  ten  larger  than  the  one-way  delay  of  the  line,  one  can  model  the  line 
as  though  it  is  electrically  short. 

For  the  model  which  we  will  develop,  we  will  assume  that  the  line  is 
electrically  short.  In  other  words,  the  time-domain  constraint  in  (3.5) 
will  place  a  bound  on  the  signals  which  may  be  analyzed  with  this  method. 


3.1.  Inductive  and  Capacitive  Coupling  Models 


board  approximate  this  situation  to  some  degree  with  the  primary  difference 
being  that  the  medium  surrounding  the  conductors  for  a  PCB  is  inhomogeneous 
(glass-epoxy,  air).  Nevertheless  this  model  provides  considerable  insight 
into  the  PCB  crosstalk  mechanism,  and  moreover,  the  resulting  transmission 
line  equations  are  solvable  in  literal  form  [19]. 

The  general,  physical  configuration  is  shown  in  Fig.  3.1.  One  conductor 
the  generator  conductor,  with  the  reference  conductor  forms  the  generator 
circuit.  This  circuit  is  driven  at  the  left  end  with  a  (pulse)  source 
Vg(t)  having  a  source  resistance  and  is  terminated  to  the  reference 
conductor  at  the  other  end  with  a  resistor  R^.  Another  conductor,  the 
receptor  conductor  with  the  reference  conductor  forms  the  receptor  circuit. 
This  circuit  is  terminated  at  the  end  near  the  source  (the  near  end)  in 
a  resistor  R^  and  at  the  end  opposite  the  source  (the  far  end)  in  a 
resistor  Rp^.  The  object  here  is,  given  the  physical  dimensions  and 
properties  of  the  conductors  and  the  surrounding  medium,  the  terminal 
resistors  and  the  time-domain  pulse  source,  Vg(t),  predict  the  time-domain 
induced  voltages  and  vF£(t). 

This  is,  of  course,  a  special  case  of  the  general,  (n+l)  conductor 

uniform  line  discussed  previously  where  n-2.  The  generator  and  receptor 

line  voltages  with  respect  to  the  reference  conductor  are  denoted  as 

v  (x,t)  and  v_(x,t),  respectively.  The  generator  and  receptor  line  currents 

are  denoted  as  i  (x,t)  and  i  (x,t),  respectively.  Again  it  is  assumed 
G  K 

that  these  currents  return  through  the  reference  conductor  such  that 
the  net  current  in  the  +x  direction  at  any  point  on  the  line  is  zero. 

The  per-unit-length  transmission-line  model  is  shown  in  Fig.  3.2. 

Once  again  the  per-unit-length  self  inductances  (capacitances)  are  denoted 
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as  i  9  £  (c^,  cD).  The  per-unit-length  mutual  inductance  (capacitance) 

G  R  G  K 

is  denoted  as  2-  (c  ).  The  transmission  line  equations  can  be  obtained 
mm  n 

from  Fig.  3.2  in  the  limit  as  ix  +  0  as 
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3i  (x,t)  9v  (x,t)  9v  (x,t) 

K  G  (  \  R 

9x  cm  9t  CR+Cm  9t 


( 3 . 6d) 


Direct  solution  of  (3.6)  in  the  time-domain  and  incorporation  of  the 
terminal  constraints  is  quite  difficult  [1].  However,  an  exact  literal 
solution  can  be  obtained  for  sinusoidal,  steady-state  excitation.  For 
sinusoidal,  steady-state  excitation,  (3.6)  becomes  a  set  of  ordinary 
differential  equations  [1]: 
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( 3 . 7d) 


where  the  phasor  line  voltages  and  currents  are  V  ,  V  ,  I  ,  I  and 

G  R  G  R 


v  v 


vG(x,t)  -  VR(x)eJu^ 


(3.8a) 


vR(x,t)  »  VR(x)e^Wt  (3.8b) 

iG(x,t)  -  IG(x)e^<iJt  (3.8c) 

iR(x,t)  -  IR(x)e^Wt  ( 3 . 8d) 

The  terminal  constraints  at  the  ends  of  the  line,  x  =  0  and  x  =  t, 

are 

V0)  -  vs  -  RsV05  (3.9a) 

V0)  “  '  rneir(0)  (3‘9b) 

VQU)  -  R^tt)  (3.9c) 

VR«)  -  R^OO  ( 3 . 9d) 


where  the  source  is  assumed  to  be  sinusoidal: 


vs(t) 


(3.10) 


The  exact  solution  of  (3.7)  with  the  terminal  constraints  in  (3.9) 
incorporated  was  obtained  in  literal  form  in  [19]  for  homogeneous  media. 
The  terminal  voltages  of  the  receptor  circuit  are 
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The  coupling  coefficient  k  is 


'('0*C»)<CR*C.) 


and  the  characteristic  impedances  of  the  respective  circuits  are 


(3.13) 


RCG  "  V  X"k 
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The  terms  V  and  I  are  the  generator  line  voltage  and  current  computed 

gdc  gdc 

from  DC  considerations: 


V  *  - =_  v 

gdc  rs+rl  s 


(3.15a) 


Igdc  "  V\  Vs 


(3.15b) 


•» :  v’-v-v  ' •>:*->££ 

'■y  -tV.v  V*V."  \  ■ 

_ *  ~*  •  '..iv  < "N-v  0.’ 


(3.19) 


Den  “  (C  ♦  ju)ST  )(C  +  jwSl  ) 
G  K 


If  the  line  is  electrically  short,  <<  X,  then 

C  =  1 
S  =  1 

The  results  in  (3.11)  simplify  to 
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If,  in  addition  to  being  electrically  short,  the  frequency  of  excitati  > 

is  small  enough  such  that  wT  «1  and  u)T  «1  then  (3.20)  simplifies  to 
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These  induced  crosstalk  voltages  can  be  calculated  from  the  equivalent 

receptor  circuit  shown  in  Fig.  3.3(a).  There  are  two  contributions  to  the 

induced  voltage.  One  source,  the  voltage  source,  depends  on  the  mutual 

inductance  between  the  two  circuits,  £  ,  and  is  referred  to  as  the  inductive 

m 

coupling  contribution.  The  other  source,  the  current  source,  depends  on 
the  mutual  capacitance  between  the  two  circuits,  c^,  and  is  referred  to  as 
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the  capacitive  coupling  contribution.  The  net  induced  crosstalk  is  the  sum 
of  the  inductive  and  capacitive  coupling  contributions. 


Both  the  inductive  and  capacitive  coupling  contributions  vary  linearly 
with  frequency  (20  dB/decade)  up  to  some  upper  frequency  limit,  f ^ .  In 
addition,  the  phase  angle  is  +90°  represented  by  j  =  /-T.  This  type  of 
behavior  is  typical  for  crosstalk.  The  upper  limit,  f  ,  is  simply  the  point 
at  which  the  increase  deviates  from  20  dB/decade.  This  was  shown  to  be 
a  function  of  the  terminal  impedances  in  [21].  Nevertheless,  this  frequency 
range  can  cover  a  large  portion  of  the  spectrum  of  the  input  pulse,  and  its 
extent  can  be  determined  either  from  frequency  response  measurements  of 
the  actual  line  (with  loads)  or  from  calculation  with  (3.11).  It  is  this 
simple,  low-frequency  region  which  we  will  concentrate  on  modeling. 


3.2.  A  Simple,  Time-Domain  Crosstalk  Prediction  Model 

One  can  determine  the  time-domain  crosstalk  pulse  using  frequency-domain 
methods  as  outlined  in  Chapter  II.  Determine  the  frequency-domain  crosstalk 
transfer  function  (magnitude  and  phase)  via  calculation  or  measurements  and 
the  spectrum  of  the  input  pulse,  v  (t),  (magnitude  and  phase)  via  Fourier 
methods.  Multiply  the  magnitudes  of  these  spectra  together  to  yield  the 
magnitude  of  the  crosstalk  pulse  spectrum,  and  add  the  phases  of  these 
spectra  to  obtain  the  phase  of  the  crosstalk  pulse  spectrum.  Then  obtain, 
via  the  inverse  Fourier  transform,  the  time-domain  crosstalk  pulse. 

This  frequency-domain  method,  although  potentially  exact,  has  a  number 
of  drawbacks.  It  is  not  suitable  for  direct  implementation  in  CAD  codes, 
cannot  handle  nonlinear  terminations,  and  requires  computation  at  a 
potentially  large  number  of  frequency  points.  It  would  be  desirable  to 
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bypass  this  time-domain  to  frequency-domain  to  time-domain  conversion 


process  and  compute  the  crosstalk  pulse  directly  in  the  time  domain  based 
on  the  time  domain  input  pulse.  This  will  be  our  objective  along  with 
any  simplification  of  that  process. 

First  let  us  suppose  that  the  frequency-domain  transfer  function 
varies  linearly  with  frequency  and  has  a  ±90°  phase  angle: 


Vf)  ■  J“  “he  Vf)  <3-22a) 

VFE(f)  -  ju  Mfe  Vs(f)  (3.22b) 


Then  it  is  a  simple  matter  to  show  that  the  time-domain  crosstalk  voltage 
is  simply  each  M  multiplied  by  the  time  derivative  of  VgCt): 

dv  ( t ) 
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A  typical  response  for  trapezoidal  Vg(t)  is  shown  in  Fig.  3.4.  The  slopes 
of  the  transfer  functions,  and  Mpg.  can  be  found  either  from  measurement 

or  from  calculation.  From  Fig.  3.3  we  obtain 
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The  near  end  and  far  end  crosstalk  pulses  will  also  be  periodic  and 
represented  by 
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Substituting  (3.22)  into  (3.28)  yields 
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where  we  have  substituted  oj=nu)o.  Comparing  (3.29)  and  (3.27)  yields  the 
result  in  (3.23). 


3.3.  An  Improved,  Time-Domain  Crosstalk  Prediction  Model 

In  the  previous  section  we  made  a  major,  simplifying  assumption. 

The  frequency  domain  transfer  function  was  assumed  to  vary  linearly  with 
frequency  and  have  a  ±90°  phase  angle  for  all  frequencies.  Clearly  this 
is  not  realistic  although,  as  we  will  show  in  the  next  chapter,  this 
very  simple  model  can  yield  quite  accurate  results  for  a  wide  range  of 
practical  problems  involving  PCB's  as  well  as  crosstalk  in  cables. 

An  important  secondary  effect  is  that  of  "common  impedance  coupling" 
as  illustrated  in  Fig.  3.5.  Consider  DC  excitation  of  the  line,  i.e., 
to»0.  The  generator  circuit  sees  a  net  resistance  of 
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Fig.  3.5.  Illustration  of  common  impedance  coupling 
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where  is  the  total  DC  resistance  of  the  reference  conductor  which  is 
assumed  to  be  much  less  than  and  RuE+^pg*  Thus  the  DC  generator  circuit 
current  is 


gdc  rs+rl 
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The  fraction  of  this  current  which  returns  through  the  reference  conductor 


is,  by  current  division. 
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Thus,  for  practical  load  resistance  values,  virtually  the  entire  generator 
circuit  current  returns  in  the  reference  conductor. 

This  results  in  a  voltage  drop  across  the  reference  conductor  of 


vo  -  Vo 


This  voltage  is  divided  across  the  receptor  circuit  loads  as 


(3.33) 
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This  induced  voltage  is  due  to  the  impedance  of  the  reference  conductor 

being  common  to  both  circuits  and  is  referred  to  as  common  impedance 

coupling.  It  tends  to  be  the  dominant  coupling  mechanism  at  the  lower 

frequencies  where  it  dominates  the  electromagnetic  coupling  computed  assuming 

perfect  conductors  as  shown  in  Fig.  3.5(b).  Depending  on  (1)  the  type  of 

reference  conductor  and  (2)  the  bit  rate  or  period  of  the  signal,  it  may 

represent  an  important  coupling  mechanism.  If  the  fundamental  (f  -1/T  ) 

o  o 

frequency  of  the  generator  signal  occurs  below  the  frequency  at  which  the 

two  coupling  mechanisms  are  equal,  common  impedance  coupling  cannot  be 

neglected.  In  the  case  of  ribbon  cables  and  data  baud  rates  below  10  kHz, 

this  was  shown  to  be  an  important  consideration  [2,3]. 

Another  important  consideration  occurs  at  the  higher  frequencies. 

The  crosstalk  does  not  continue  to  increase  at  a  20  dB/decade  rate  but 

replicates  causing  nulls  and  maxima  to  appear  in  the  frequency  response.  The 

magnitude  of  the  maxima  and  the  locations  of  the  nulls  are  difficult  to 
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predict  but  occur  above  the  point  at  which  the  line  becomes  4  X  in  length. 
Other,  secondary  minima  (not  nulls)  can  occur  below  this  point  [20].  We 
propose  to  bound  this  high  frequency  region  with  a  constant  coupling  level 
rather  than  have  it  increase  with  frequency  indefinitely.  Thus,  above  f^  in 
Fig.  3.3  we  choose  a  constant  coupling  level  at  the  level  existing  at  f  and 
zero  phase  angle. 

Combining  these  two  practical  considerations  yields  the  proposed 
frequency-domain  crosstalk  model  shown  in  Fig.  3.6.  The  curve  is  defined  by 
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Substituting 
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(3.41) 


yields 
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(3.42) 


For  a  general  v  (t)  this  requires  direct  integration  (or  convolution). 
However  if  that  is,  the  frequency  response  is  assumed  to  increase  at 

20  dB/decade  indefinitely  but  common  impedance  coupling  is  retained,  we 
obtain 


or 


-  K,  <1  .  -f> 

yS(t>  ^  "l 


v(t)  -  vg(t)  ♦  M 


dvs(t) 

dt 


(3.43) 


(3.44) 


as  in  the  previous  section. 


An  alternative  derivation  of  (3.42)  can  be  obtained  by  viewing  the 
frequency  response  given  in  Fig.  3.6  as  an  analogue  to  that  obtained  by  a 
high-pass  filter.  This  is  not  to  imply  a  physical  correspondence  but  only 
a  mathematical  analogue  for  the  purposes  of  deriving  (3.42).  For  example, 
the  high-pass  filter  shown  in  Fig.  3.7  has  the  frequency  response  transfer 
funct ion 
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Comparing  (3.45)  and  (3.36)  shows  that,  for  this  analogue, 
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(3.45) 
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(3.46b) 


Wu  cTr“|  |  (Ra+Rb)] 


(3.46c) 


Thus,  if  one  desires,  the  circuit  shown  in  Fig.  3.7  may  be  used  to  simulate 
the  t ime-domain  line  response  once  the  analogous  quantities  R,  R^,  R^,  C  are 
calculated  from  measured  (or  calculated  [19])  results  from  the  following. 
There  are  three  independent  parameters  in  Fig.  3.6,  for  example,  K^,  u)^ 
and  M.  The  circuit  in  Fig.  3.7  contains  four  parameters,  R,  C,  R^,  R^. 

So  we  can  constrain  one  and  find  the  remaining  three.  A  number  of 
possibilities  exist.  One  may  choose  C  arbitrarily  and  compute  R,  R^  and  R^ 
from  (3.46).  If  we  do  so  we  find,  using  M  -  K^/w^ 
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S  V  ’ 


Fig. 


.7.  A  lumped,  equivalent  circuit  for  simulation  of 
time-domain  crosstalk. 
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Note  that  K  =  Mo  .  Since  K  >  K,  the  denominators  of  R  and  R^  are  positive 
u  u  u  L  a  d 

However  it  is  possible  to  have  >  1.  This  would  require  R^  to  be  negative. 
Most  CAD  codes  can  handle  negative  resistors. 

Thus  one  can  use  the  lumped  circuit  of  Fig. 3. 7  in  any  standard,  lumped- 
circuit  CAD  code  to  calculate  the  induced,  time-domain  crosstalk  pulses  once 
the  circuit  parameters  are  calculated  with  (3.47)  from  frequency-domain 
measurements.  It  should  be  stressed  again  that  no  physical  correspondence 
between  the  circuit  in  Fig.  3.7  and  the  actual  transmission  line  are  intended 
nor  should  be  attempted.  The  circuit  of  Fig.  3.7  is  only  intended  as  an 
analogue  which  provides  a  convenient  way  of  calculating  the  line  time-domain 
crosstalk  using  a  CAD  code. 

For  known  v  (t),  use  of  the  equivalent  circuit  is  not  necessary.  In  the 
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next  section  we  will  give  an  explicit  solution  to  (3.42)  for  a  specific 
v  (t) — a  periodic,  trapezoidal  pulse  train — which  is  common  on  PCB's. 

O 

3.3.1.  Solution  for  Trapezoidal  (Clock)  Pulses 

Consider  a  periodic,  trapezoidal  pulse  train  shown  in  Fig.  1.1(a).  This 
type  of  signal  is  commonly  present  on  PCB  lands  and  is  meant  to  resemble 
clock  signals.  For  this  type  of  pulse  we  may  obtain  an  explicit  solution 
to  (3.42).  We  will  assume  that  the  line  is  initially  uncharged  when  the 
pulse  is  applied  at  t»0,  v(0)=0.  Thus  (3.42)  becomes 
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The  result  requires  a  determination  of  the  convolution  of  e  UJut  and  two 
signals,  vg(t)  and  vg(t): 
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(3.49a) 


(3.49b) 


where  *  denotes  convolution.  The  signal,  v  (t),  and  its  derivative  v  (t), 

b  S 

are  shown  in  Fig.  3.8.  The  convolutions  in  (3.49)  are  shown  in  Fig.  3.9. 

Each  convolution  requires  the  product  of  e  and  either  v  (X)  or  v  (X) 

O  b 

as  t  progresses  and  the  determination  of  the  resulting  area  under  the  product 
curve.  This  interpretation  shows  three  distinct  regions  over  a  complete 
period: 
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and  x  is  the  pulse  width  (between  50%  points). 


One  of  the  simplest  ways  of  determining  the  solution  to  (3.48)  for  this 

signal  is  to  recognize  that  v_(t)  is  a  sequence  of  pulses  and  may  be  written 

5 

in  terms  of  the  unit  step  function  u(t)  over  0  -  t  -  T  as 

vs(t)  -  J-  [u(t)  -  u(t-tl)] 
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—  [u(t-t_)  -  u(t-t-)] 
Tf  J 


(3.52) 
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Similarly,  v,,(t)  can  be  written  in  terms  of  unit  ramp  functions,  r(t),  over 
0  -  t  -  T  as 


vg(t)  =  [r(t)  -  r(t-t^) ] 

*J  r 


-  qr-  [r(t-t_)  -  r(t-t-)] 
Tf 
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(3.53) 


We  will  determine  the  basic  convolutions  in  (3.49)  in  response  to  a  unit 
step  or  a  unit  ramp  and  use  (3.52)  and  (3.53)  to  determine  the  complete 
response . 

The  convolution  of  e  Wut  and  a  unit  step  is 
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The  convolution  of  e  U)ut  and  a  unit  ramp  is 
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(3.55) 


Note  that  (3.55)  is  the  integral  of  (3.54)  as  it  should  be. 

Thus,  the  solution  to  (3.48)  is  the  appropriately  weighted  combination 
of  these  solutions  delayed  according  to  (3.52)  and  (3.53): 
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Although  these  equations  appear  to  be  complicated  they  may  be  broken 
down  according  to  the  appropriate  time  intervals  of  interest. 


For  0  <  t  $  t^  ,  (3.56)  becomes 
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0  <  t  <  t 


The  first  term  is  simply  the  input  signal  multiplied  by  the  common  impedance 
coupling  factor,  IC  ,  which  is  sensible  to  expect.  Since  K  -  Mm  and 
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>>  for  typical  boards,  the  second  term  may  be  simplified  so  that 


r(t)  -  KjV^t)  +  —  (1  -  e~<i,ut:) 


(3.58) 


The  second  term  is  the  result  derived  in  the  previous  section  where  we 
neglected  common  impedance  coupling  and  assumed  the  frequency  response  to 
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increase  indefinitely  with  frequency  but  is  multiplied  by  the  term  (1  -  e  u 
For  1  t  1  tj  additional  terms  are  added.  The  result  is 
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But  this  may  be  written  as 
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Since  we  usually  have  >>  K^,  this  simplifies  to 

v(t)  -  VKL  +  ^  1 1  -  e"U,uTr]  «-Mt“Tr> 
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The  first  term  is  the  input  signal  over  this  interval  multiplied  by  the 
common  impedance  coupling  coefficient.  The  second  term  is  the  initial 
condition  at  t  ■  tj  ■  resulting  from  the  second  term  of  (3.58)  multiplied 
by  the  decay  term  e  Tr^: 
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this  can  be  simplified  to 
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For  £  t  £  T  the  result  is 
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The  first  term  is  part  of  the  result  due  to  the  initial  condition  established 
on  the  previous  interval  as  is  the  third  term  which  decays  by  e-U>u^t  . 
Since  >>  we  may  again  replace  (M  -  with  M.  Note  that  the  usual 

term  involving  the  product  of  and  v,,(t)  over  this  time  interval  is  absent 
since  v<,(t)  ■  0  over  this  time  interval. 

In  summary,  the  response  consists  of  three  contributions: 


v(t)  ■  K^v,, ( t )  +  Mf^(t)  +  Mf2(t) 


(3.67a) 


The  function  f ^ ( t )  is  the  response  of  a  first-order  circuit  with  time  constant 
l/u>u  to  a  square  pulse  of  height  v/xr  occuring  over  0  5  t  S  t^.  The  function 
f 2 ( t )  is  the  response  of  a  first-order  circuit  with  time  constant  l/<i>u  to  a 
square  pulse  of  height  -  V/x^  occurring  over  tj  1  t  $  tj.  Thus 
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3.4.  Implementation  in  CAD  Codes 

It  should  be  emphasized  that  the  intent  in  developing  a  model  for  use  in 
CAD  codes  is  to  obtain  one  which  (1)  provides  estimates  of  the  time-domain 
crosstalk  and  (2)  does  not  consume  an  excessive  amount  of  computation  time 
or  array  storage.  The  complexity  of  modern  PCB's  provides  strong  motivation 
for  the  latter  criterion. 

The  simplest  model  would  be  the  one  addressed  in  section  3.2  in  which 

we  assumed  -*•  i.e.,  the  frequency  response  of  the  transfer  function 

increased  linearly  with  frequency  indefinitely.  To  implement  this  model  we 

note  that  ju)  in  Fig.  3.3  may  be  replaced  by  and  the  result  in  (3.23) 

will  again  be  obtained.  Thus  equivalent  sources  may  be  inserted  into  all 

receptor  circuits  and  the  values  of  these  equivalent  sources  depend  only 

on  land  geometry  (through  l  and  c  )  and  the  derivatives  of  the  signals 

mm 

on  the  other  lands.  This  is  illustrated  in  Fig.  3.12.  In  this  way,  nonlinear 
load  terminations  may  be  handled  since  we  only  model  the  induced  sources  in 
the  receptor  line. 

The  degree  to  which  this  approximation  is  valid  is  usually  difficult 
to  determine.  Nevertheless  it  represents  a  simple  approach  to  estimating 
the  solution  to  a  difficult  problem. 

The  question  now  arises  as  to  whether  it  is  worth  implementing  the 
improved  model  of  section  3.3.  In  the  author's  opinion  the  answer  is  no. 

To  do  so  would  require  multiplication  of  each  term  in  each  source  by  expon¬ 
ential  decay  terms.  Each  would  depend  on  the  generator  circuit  and  would 
probably  require  more  preliminary  analysis  than  analysts  would  be  willing 
to  expend.  In  addition  each  oj^  depends,  in  addition  to  readily  available 
parameters  such  as  land  geometry  (through  l  and  c  ) ,  on  the  loads  of  the 
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receptor  circuit-  Thus  one  could  not  simply  insert  sources  which  depended 
only  on  land  geometry. 

Also  it  should  be  kept  in  mind  that  the  simple  model  in  Fig.  3.12 
provides  an  upper  bound  on  the  crosstalk.  Generally,  one  is  interested  in 
CAD  codes  in  obtaining  estimates  of  potent ial  crosstalk  problems.  Culling 
the  large  number  of  those  cases  which  will  not  provide  crosstalk  problems 
will  be  a  significant  contribution  since  the  analyst  can  devote  additional 
effort  to  determining  whether  those  suspect  cases  in  fact  would  provide 
problems . 

When  all  of  these  factors  are  considered,  it  seems  that  the  simple 
model  in  Fig.  3.12  provides  the  logical  choice  for  CAD  implementation. 


3.5.  A  Further  Improvement 

The  models  of  the  previous  sections  represent  approximations  to  the 
exact  solution  given  in  [19]  and  outlined  in  section  3.1.  It  is  possible 
to  find  situations  which  do  not  satisfy  the  restrictions  of  those  models 
even  though  the  line  is  electrically  short.  For  example,  there  exist  cases 
which,  even  for  resistive  loads,  the  frequency  response  does  not  vary  at 
20  dB/decade  in  the  region  between  f  and  f  in  Fig.  3.6  [20,21].  In  order 
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to  handle  these  cases,  we  use  the  development  given  in  [20]. 

The  exact  solution  in  [19]  can  be  written  as  [20] 
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and  and  are  given,  as  before,  by  (3.24).  The  other  terms  in  (3.69) 
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where  T  and  T  are  given  by  (3.18)  and 
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Suppose  that  we  only  consider  frequencies  of  the  excitation  pulse. 


vc(t),  such  that  the  line  is  electrically  short,  X.  Then  we 
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Replacing  jto  with  the  differential  operator  p  we  obtain  the  differential 

equations  relating  v  (t)  and  v  (t)  to  the  input  voltage  v  (t)  as 
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(3.74a) 
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These  differential  equations  can  be  solved  for  v  (t)  and  v  (t)  for 
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t  >  0  given  the  initial  conditions  at  t=0: 
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(3.75) 


This  form  of  solution  can  be  used  assuming  that  the  initial  conditions 

in  (3.75)  are  all  zero.  Perhaps  the  assumption  v„,„(0)  »  v„„(0)  -  0  is 

-  NE  FE 

reasonable,  but  assuming  the  derivatives  to  be  zero  at  t=*0  is  not  an  obviously 
valid  result. 

Nevertheless,  Lf  we  assume  all  initial  conditions  are  zero  and  write 
the  equations  in  terms  of  their  Laplace  transforms,  we  obtain 
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If  the  two  lines  are  weakly  coupled,  \i>  <<  1,  then 
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Thus  the  denominators  of  (3.76)  factor  into  a  product  of  two  factors  each 
of  which  involve  the  time  constant  of  one  circuit.  Suppose  we  assume  that 
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The  Laplace  transform  of  this  is 
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The  transfer  functions  in  (3.76)  become 
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Vtr  -  V  <•  *  1/tr> 


VFE(S) 


i 

:tg  -  V  (s  *  u\ 


<tr  -  tg)  (s  *  1/tr> 


The  inverse  Laplace  transforms  become 

,  .  r  (t '  V'"t/Tc  <T  -  Ve‘t/T“  V 

Vne  tg  '  tr  (,r  -  tg> 

»_(t)  -  t  -e-t/T0  .  T 


(3.78) 


(3.79) 


(3.80a) 


(3.80b) 


(3.81a) 


(3.81b) 


(3.82a) 


(3.82b) 


The  solution  is  therefore  the  familiar  double  exponential  pulse  similar  to 
the  overdamped  solution  of  a  second  order  circuit. 

For  practical,  trapezoidal  pulses  typical  of  clock  pulses,  the  solution 
is  more  complicated  than  for  methods  of  the  previous  sections.  Thus  this 
additional  refinement  is  not  suitable  for  CAD  implementation  unless  (1)  a 
suitable  analog  circuit  can  be  derived  and  (2)  the  initial  conditions  in 
(3.75)  are  known  and  translated  to  appropriate  initial  conditions  on  this 
analog  circuit. 


IV.  Experimental  Results 

Experimental  results  were  obtained  to  illustrate  the  prediction 
accuracies  as  well  as  the  restrictions  of  the  simple  models  developed  in  the 
previous  chapter.  Two  types  of  lines  were  tested — a  typical  PCB  and  a 
ribbon  cable.  Although  the  primary  intent  of  this  report  is  to  address 
crosstalk  on  PCB's,  it  is  worthwhile  to  show  that  these  models  can  predict 
time-domain  crosstalk  in  fairly  long  cables  so  long  as  the  pulse  rise/fall 
times  are  sufficiently  long  as  discussed  previously. 

This  chapter  will  be  organized  into  a  discussion  of  (1)  the  line 
characteristics,  (2)  the  frequency  response  data,  and  (3)  the  time-domain 
response  data.  Model  predictions  for  both  frequency-domain  and  time-domain 
results  will  be  shown  using  (1)  the  lumped-circuit  iterative  models  and  (2) 
the  simple  prediction  model  of  the  previous  chapter. 

The  lines  will  be  driven  and  terminated  as  in  Fig.  3.1.  The  input 
signal  v  (t)  will  be  monitored  and  set  at  the  desired  level.  Thus  R  =0. 

The  other  three  terminal  resistors  will  be  chosen  equal: 

^  “  RNE  “  rfe  =  R 
Two  values  of  R  will  be  used 

R  -  500,  1  kQ 

These  two  values  were  chosen  for  the  investigation  in  order  to  enhance  either 
inductive  or  capacitive  coupling.  We  will  find  that  both  the  frequency- 
domain  and  the  time-domain  crosstalk  properties  depend  quite  strongly  on 
whether  inductive  or  capacitive  coupling  is  dominant. 
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The  measurement  and  excitation  equipment  are  as  follows.  For  the 
frequency-domain  results,  the  sources  were  standard  sinusoidal  oscillators 
which  covered  the  measurement  range  of  10  kHz  to  500  MHz  for  the  PCB  and 
100  Hz  to  100  MHz  for  the  ribbon  cable.  The  measurement  equipment  for  the 
frequency-domain  results  consisted  of  HP  3400A  RMS  voltmeters  for  use  below 
1  MHz  and  a  HP  8405A  vector  voltmeter  for  use  above  1  MHz. 

For  the  time-domain  measurements,  the  excitation  source  was  a  HP  8015A 
pulse  generator.  The  crosstalk  voltages  were  measured  with  a  Tektronix  7834 
oscilloscope.  For  the  PCB,  P6201  active  FET  probes  were  used.  Several 
values  of  rise/fall  time  were  used.  For  the  PCB,  a  1  MHz  trapezoidal  wave¬ 
form  which  transitioned  from  -1.25  V  to  +  1.25  V  was  used.  The  rise/fall 
times  of  the  pulses  which  were  investigated  were  10  ns,  50  ns  and  100  ns. 

The  board  one-way  delay  was  on  the  order  of  1.3  ns.  Thus  we  will  investigate 
rise/fall  times  which  approach  the  limit  discussed  in  Chapter  2. 

For  the  ribbon  cable,  a  20  kHz  trapezoidal  waveform  typical  of  RS-232C 
applications  was  used.  The  pulse  transitioned  from  -1.25  V  to  +  1.25  V 
and  had  rise/fall  times  of  6  ps,  1  ps,  700  ns,  400  ns,  100  ns,  60  ns  and 
20  ns.  The  line  one-way  delay  was  approximately  32  ns. 

4.1.  Printed  Circuit  Board 

The  printed  circuit  board  consisted  of  glass-epoxy  substrate  62  mils  in 
thickness  with  a  ground  plane  on  one  side  and  lands  on  the  other.  A  photo¬ 
graph  is  shown  in  Fig.  4.1.  All  conductors  are  1  ounce  copper  (1.4  mils 
thickness).  The  lands  are  100  mils  in  width  and  separated  by  100  mils.  The 
land  lengths  are  20  cm  Q£-20  cm) . 


4.1.1.  Frequency-Domain  Results 


The  frequency-domain  crosstalk  was  measured  from  10  kHz  to  500  MHz. 

The  results  are  shown  for  near  end  crosstalk  and  R-50Q  in  Fig.  4.3,  for  far 
end  crosstalk  and  R-50Q  in  Fig.  4.4,  for  near  end  crosstalk  and  R»i  kQ  in 
Fig.  4.5,  and  for  far  end  crosstalk  and  R**l  kQ  in  Fig.  4.6.  The  predictions 
of  the  multiconductor  transmission  line  model  are  also  shown  [1,2]. 

From  these  results  we  obtain 


Table  4.1 


R 

M 

0) 

u 

near 

(  50Q 

0 

1.06E-10 

2E9 

end 

1  1  kQ 

0 

6.37E-10 

6.28E7 

far 

(  500 

0 

-4.46E-11 

7E9 

end 

1  1  kQ 

0 

6.37E-10 

6.28E7 

The  predictions  of  the  lumped  Pi  and  Tee  models  for  1  and  5  sections 
are  shown  in  Fig.  4.7  through  Fig.  4.14.  These  are  compared  against  the 
predictions  of  the  multiconductor  transmission  line  model.  Generally  the 
lumped-circuit  iterative  models  provide  accurate  predictions  up  to  the  point 

_L 

at  which  the  line  is  10  A  (approximately  80  mHz).  Note  that  increasing  the 
number  of  lumped-circuit  iterative  model  sections  does  not  result  in  a 
dramatic  increase  in  frequency  coverage. 


'.'i 


4.1.2.  Time-Domain  Results 

A  1  MHz  trapezoidal  waveform  which  transitioned  from  -1.25  V  to  1.25  V 
was  applied.  Three  rise/fall  times  were  investigated — 100  ns,  50  ns,  10  ns. 
The  board  one-way  delay  was  1.3  ns.  Thus  the  10  ns  rise/fall  time  represent! 
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Fig.  4.5.  Frequency  response  experimental  data  vs. 

predictions  of  the  MTL  model.  Near  End,  1  kQ. 
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Fig.  4.6.  Frequency  response  experimental  data  vs. 

predictions  of  the  MTL  model.  Far  End,  1  kQ 
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Lumped  iterative  model  vs.  MTL  model  predictions. 
Near  End,  Pi,  1  kQ. 


osstalk,  ZL=  1  K 


Table  4.3 


Maximum  Crosstalk  Levels 
(R-l  kO) 


Predicted 

Measured 

VNE 

12.7  mV 

14  mV 

T  -T,-100  ns 

VFE 

12.7  mV 

14  mV 

VNE 

24.4  mV 

25  mV 

T  «t,«50  ns 

r  i. 

VFE 

24.4  mV 

25  mV 

VNE 

59.4  mV 

55  mV 

T  -T  -10  ns 

L  L 

V 

FE 

59.4  mV 

65  mV 

The  predicted  maximum  levels  are  remarkably  accurate.  The  exponential 

|P 

factor  (1-e  )  was  of  no  consequence  for  R-50Q.  The  experimental  results 

show  this — the  crosstalk  pulses  are  virtually  rectangular  and  show  no  signi¬ 
ficant  rounding  of  the  corners.  Thus  for  R-50Q  one  can  predict  the  cross¬ 
talk  as  simply  the  slope  of  the  frequency  response,  M,  multiplied  by  the 
slew  rate  of  Vg(t).  However,  this  factor  was  important  for  all  of  the 
R-l  kQ  results.  The  experimental  results  show  this  in  that  all  pulses 
show  exponential  rise  and  decay. 

Note  that  all  crosstalk  pulses  have  the  same  polarity  as  the  slope  of 
the  input  signal  except  for  R-50Q  and  far  end  crosstalk.  It  was  pointed 


out  in  Chapter  3  that  this  would  be  the  case  if  inductive  coupling  dominated 


capacitive  coupling  so  that  M^g  is  negative. 


The  predictions  of  the  lumped-circuit  iterative  Pi  and  Tee  models  using 


5  sections  are  shown  in  Fig.  4.21  through  Fig.  4.32.  The  predictions  of 


the  simple  model  using  the  frequency-domain  parameters  of  Table  4.1  and  the 


circuit  analogue  in  Fig.  3.7  are  also  shown  on  these  figures.  The  simple 


model  provides  reasonable  predictions  and  requires  considerably  less 


computational  expense  than  the  lumped  iterative  models. 


4.2.  Ribbon  Cables 


A  three-wire  ribbon  cable  was  tested.  All  wires  were  #28  AWG  (7x32) 


and  the  line  length  was  4.737  m  (15'6V).  The  dielectric  insulation  was 


polyvinyl  chloride.  The  wires  are  separated  (center-to-center)  by  50  mils. 


The  center  wire  was  chosen  as  the  reference  conductor  and  the  other  wires 


were  terminated  as  shown  in  Fig.  4.33.  Once  again  Rg-0  and 


*L  "  ‘  R 


with  two  values  of  R  being  used,  R-50Q,  1  kQ. 


4.2.1.  Frequency  Domain  Results 


The  frequency  domain  crosstalk  was  measured  from  100  Hz  to  100  MHz. 


The  results  are  shown  for  near  end  crosstalk  and  R*50Q  in  Fig.  4.34,  for  far 


end  crosstalk  and  R=50Q  in  Fig.  4.35,  for  near  end  crosstalk  and  R=1  kQ 


in  Fig.  4.36,  and  for  far  end  crosstalk  and  R-l  kQ  in  Fig.  4.37.  The 


predictions  of  the  multiconductor  transmission  line  model  are  also  shown 


[1,2,3]. 
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Fig.  4.26.  Time-domain  predictions 
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Fig.  4.27.  Time-domain  predictions 
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4.29.  Time-domain  predictions 

1  kQ,  Near  End,  T  “  Tf=  50  ns 
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Fig.  A. 34.  Frequency  response  experimental  data  vs. 

predictions  of  the  MTL  model.  Near  End,  500 
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From  these  results  we  obtain 


Table  4.4 


8.9E-3 

1.26E-8 

4 . 5E-4 

2.16E-8 

8.9E-3 

-1  E-8 

4 . 5E-4 

1.68E-8 

The  predictions  of  the  lumped-circuit  iterative  Pi  and  Tee  models  are 
shown  in  Fig.  4.38  through  4.45.  These  models  provide  accurate  predictions 

_L 

for  one  section  up  to  a  frequency  where  the  line  is  on  the  order  of  10  X 
in  length  (32  MHz). 

4.2.2.  Time-Domain  Results 

A  20  kHz  trapezoidal  pulse  train  typical  of  RS-232C  data  transmission 
was  applied  to  the  cable.  The  pulse  transitioned  from  -1.25  V  to  +1.25  V. 
The  rise/ fall  times  used  were  6  ps,  1  ps,  700  ns,  400  ns,  100  ns,  60  ns 
and  20  ns.  The  one-way  line  delay  is  32  ns.  Thus  the  100  ns  rise/fall  time 
is  approaching  the  limit  of  applicability  of  the  model.  The  results  are 
shown  in  Fig.  4.46  through  Fig.  4.59. 

From  the  frequency-domain  results  for  M  and  w^,  we  compute  the  maximum 
crosstalk  using  the  model  in  section  3.3  as  seen  in  Table  4.5  and  Table  4.6. 

The  predicted  maximum  levels  are  quite  accurate  except  for  the  very 
short  rise/ fall  times.  For  T^-T^-10  ns  and  20  ns  we  have  more  errors  than 
encountered  with  the  PCB.  Since  the  line  one-way  delay  time  is  on  the 
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Fig.  4.42.  Lumped  iterative  model  vs.  MTL  model  predictions 
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Table  4.5 


Maximum  Crosstalk  Levels 
(R-50Q) 


Predicted 

Measured 

VNE 

4 . 2  mV 

9  mV 

T  -b  us 

r  f 

FE 

-  3 . 3  mV 

9  mV 

VNE 

25.2  mV 

28  mV 

T  «T  -I  US 
r  f 

VFE 

-  20  mV 

-  24  mV 

36  mV 

34  mV 

T  =T  *700  ns 

VFE 

-  28.6  mV 

-  30  mV 

v„„ 

63  mV 

55  mV 

T  «t  =400  ns 

VFE 

-  50  mV 

-  45  mV 

T  -T  -100  ns 

VNE 

217.9  mV 

140  mV 

VFE 

-172.9  mV 

-  110  mV 

T  -T  =60  ns 

VNE 

293.5  mV 

170  mV 

VFE 

-232.9  mV 

-  140  mV 

VNE 

415.4  mV 

200  mV 

T  =T  =20  ns 

VFE 

-329.7  mV 

-  180  mV 

order  of  32  ns,  this  is  not  surprising.  For  rise/fall  times  of  400  ns  and 

longer,  the  predictions  are  very  good.  For  400  ns  rise/fall  times  1^=12. 5 

Note  again  as  with  PCB's,  the  near  end  and  far  end  crosstalk  polarities 
are  the  same  as  the  appropriate  slope  of  the  input  pulse  except  for  R*=50Q 
and  far  end  crosstalk.  Once  again  this  is  due  to  inductive  coupling 
dominating  capacitive  coupling  so  that  is  negative. 

The  predictions  of  the  lumped-circuit  iterative  Pi  and  Tee  models  for 
5  sections  are  shown  in  Fig.  4.60  through  4.71  for  rise/fall  times  of  700  ns, 

100  ns,  and  20  ns.  The  predictions  are  reasonably  accurate.  Also  we  have 

shown  the  predictions  of  the  simple  model  using  the  measured,  frequency- 
domain  parameters  of  Table  4.4  and  the  lumped  circuit  analogue  of  Fig.  3.7. 

4 . 3  Summary 

The  lumped-circuit  iterative  models  predicted  frequency-domain  results 

_L 

for  frequencies  up  to  the  point  at  which  the  line  was  10  X  in  length.  Above 
that  point  more  sections  are  needed  but  very  little  extension  of  the  valid 
frequency  range  is  obtained. 

The  simple  model  of  Chapter  3  provided  quite  accurate  predictions  of 
the  maximum  crosstalk  levels  so  long  as  the  pulse  rise/fall  time  was 
longer  than  10  to  20  times  the  one-way  line  delay.  In  addition,  considerable 
insight  into  the  coupling  phenomenon  is  gained  with  the  simple  model.  The 
lumped-circuit  iterative  models  provide  no  such  insight. 
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The  purposes  of  this  report  have  been  (1)  to  investigate  the  ability 
of  available,  time-domain,  crosstalk  prediction  models  and  (2)  to  determine 
a  simple  time-domain  crosstalk  prediction  model  which  is  suitable  for 
inclusion  in  CAD  codes.  The  primary  type  of  transmission  line  structure 
which  these  models  were  to  address  was  the  typical  printed  circuit  board 
with  signals  which  represent  typical  bit  rates  and  rise/fall  times. 

Within  these  objectives  and  restrictions,  a  simple  model  suitable 
for  inclusion  in  CAD  codes  was  developed.  The  model  requires  two  groups 
of  information.  The  first  set  is  related  to  the  cross-sectional  structure 
of  the  lines  ar/d  are  contained  in  the  per-unit-length  mutual  inductances 
and  capacitances  between  all  lines.  The  other  set  are  the  voltages  and 
currents  on  the  interfering  lines.  The  first  set  of  data,  the  mutual 
inductances  and  capacitances  can  be  obtained  from  direct  calculation,  direct 
measurement  or  inferred  from  frequency-domain  crosstalk  measurements. 

The  second  set  of  data — the  voltages  and  currents  on  the  interfering  lines — 
are  normally  calculated  by  CAD  codes  and  can  be  used  directly.  The 
crosstalk  sources  depend  on  these  mutual  elements  and  the  slew  rates  of 
the  voltages  and  currents  on  the  interfering  lines. 

Experimental  results  for  a  typical  FCB  as  well  as  a  ribbon  cable  show 
that  the  simple  model  is  capable  of  yielding  prediction  accuracies  better 
than  3  dB  so  long  as  the  pulse  rise/fall  time  is  at  least  a  factor  of  10 
to  20  larger  than  the  one-way  line  delay.  Even  if  this  restriction  is 
exceeded  the  simple  model  provides  an  upper  bound  on  the  maximum  crosstalk. 

When  one  considers  the  complexity  of  even  relatively  simple  PCB's, 
it  becomes  quite  evident  the  sufficient  computing  power  does  not  exist 


for  analyzing  time-domain  crosstalk  between  all  lands  in  real  time  with 
sophisticated  prediction  models  or  even  the  relatively  simple  lumped-circuit 
iterative  models  (using  only  one  section  per  line).  The  simple  model 
proposed  here  has  the  potential  for  handling  this  problem. 

In  addition,  perhaps  the  minimum  amount  of  information  is  required 
to  be  gathered  by  the  analyst.  However,  it  remains  to  provide  an  inter¬ 
mediate  code  to  transform  the  physical  land  layout  data  to  the  mutual 
inductances  and  capacitances  required  by  the  model. 
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